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(1) 


Wirminghoff  [6],  in  his  doctoral  thesis  investigated  the  effect  on 
the  process  of  geostrophic  adjustment  [2]  in  one  dimension  of  several 
different  spatial  discretization  schemes.   His  results  were  later  in- 
corporated by  Arakawa  and  Lamb  [1] .   The  model  he  chose,  which  had 
been  previously  investigated  by  Rossby  [4],  Cahn  [3]  and  others,  was 
the  linearized  one-dimensional  shallow  water  equations  with  no  mean 
flow,  in  an  infinite  region: 

_  _  fv  +  g  _  „  o  (1) 

3v 

f^  +  fu  =  0  (2) 

f+Hf=0,  (3) 

where  u  is  the  perturbation  velocity  in  the  x  direction,  v   is  the 
perturbation  velocity  normal  to  the  x  direction,  H  and  h   are  the  mean 
and  perturbed  heights  of  the  free  surface,  respectively,  and   g  >  0   and 
f  >  0  are  gravitational  and  Coriolis  parameters,  respectively.   The  purpose 
of  this  note  is  to  point  out  certain  errors  in  Winninghof f ' s  solutions, 
present  the  correct  solutions,  and  reexamine  his  conclusions  in  light  of 
the  correct  solutions. 

Winninghof f  followed  earlier  studies  in  solving  (l)-(3)  by  eliminating 

between  the  equations  to  arrive  at: 

2  2 

3  u    c2  3  u  , .  v 

— -  +  f  u  -  gH  — -  =  0  ,  (4) 

3tZ  3xZ 

then  solving  (4)  by  a  Fourier  Transform  approach.   After  solving  (4), 
solutions  for  h  and  v  are  obtained  by  substitution  into  (2)  and  (3) , 


(3) 


although  closed-form  solutions  are  not  presented  in  some  of  the  papers. 
Note  the  dispersive  character  of  (4)  is  clearly  seen  by  assuming  a  wave 
solution: 


.        .    .   i(kx-vt) 
u(x,t)  =  A  e 


which  implies: 


v2  =  f2  +  gHk2  =  f2  (1  +  A2k2),  \   =  /gUTf"  .  (5) 


In  [5],  we  presented  a  method  where  (l)-(3)  were  transformed  directly, 
then  the  transformed  problem  was  solved  as  a  system  of  coupled  ordinary 
differential  equations.   If  we  denote  transforms  by  an  overhead  tilde,  e.g, 

on 

-ikx 


u(k,t)  =   /   u(x,t)e     dx  , 


then  (l)-(3)  reduce  to 


du  ~  ~ 
—  -  f v  -  lkg  h 
dt 


dt"  =  "fu  (6) 


dh  ~ 
-j—  =  -lkH  u 
dt 


Now,  if  we  denote  initial  values  by  a  subscript,  e.g, 

00 

-ikx 


u  (k)  =   /   u(x,0)e  1KX  dx  =  u(k,0)  , 


0 

—  00 

and  use  Winninghof f 's  initial  conditions: 

u(x,0)  =  0 
h(x,0)  =  0 

-a  <  x  <  a 


,0,  -  jV°'  • 

(0,  ol 


v(x, 

)therwise  , 

then  using  equation  (12)  of  [5],  the  transformed  solution  is 


(4) 


(7) 


/i  ^\       0   .  ,      s  sxn(ak) 
u(k,t)  =  sin(vt)  J — - 

!2      2        1 
k  gH    f      ,   J  sin(ak)  ,ON 

— |-  +  -^  cos(vt)) g — t.  (8) 

2ikHfV_  .  /  , x 

7"/,  -n     O/i       /  „\i  sxn(ak) 

h(k,t)  = ^ tl  -  cos(vt)}  g — -   . 

v 
Winninghoff,  using  elimination,  solved  only  for  u(x,t)  , 


(9) 


u(x,t)  =  27  /  u(k,t)e1  X  dk 

—  00 

2faV    °° 
=  0  J      sin(ak)  sin(vt)  cqs  kx  dk 

it    jL    ak       v 

then  solved  for  h(x,t)   using  the  continuity  equation.   [If  this  is  done 

to  (9),  or  if  the  third  equation  in  (8)  is  inverted  directly,  we  find 

2HfV0   - 
h(x,t)  = /  — cos^vtM  sln(ak)  sin(kx)dk  .         (io)] 

77    0       v 
The  error  in  his  solution  occurs  in  converting  from  the  differential  case 
to  the  discretized  case,  using  any  one  of  the  difference  schemes  he  denotes 
as  schemes  A-E.   (Figure  1)   He  performs  this  conversion  by  simply  replacing 
v   in  (9)  by  its  equivalent  finite  difference  form,  evaluating  the  integral 
numerically,  then  finding  h(x,t)   from  the  continuity  equation.   We  now 
show  this  approach  neglects  the  effects  of  spatial  averaging  of   u(x,t) 
and   v(x,t)   that  occurs  in  Schemes  C  and  D.   For  example,  in  Scheme  C, 
the  finite  difference  form  of  (l)-(3)  is: 

3u  m   f  [v(x+d/2,t)  +  v(x-d/2,t)]     [h(x+d/2,t)  -  h(x-d/2,t)] 
3t  2  8  d 

ll  =  _f  [u(x+d/2,t)  +  u(x-d/2,t)] 
3t  '  2 

3h    „  [u(x+d/2,t)  -  u(x-d/2,t)] 
—  -  -H  -  . 


(5) 


If  we  transform  this  system,  distinguishing  from  the  differential 
system  by  using  a  subscript   c  on  the  variables,  we  have: 

c        ^           /(/1mu          .      sin(kd/2)    ~ 
=   f   cos(kd/2)v     -   lg  ,,,,, h 


dt  *    —  v-'^vc        -s        (d/2)        »c 

dv 

-^  =  -f   cos(kd/2)uc  (12) 

c  #TT   sin(kd/2)    ~ 

=  -iH  ,  \  ,  „\ u 


dt  (d/2)  c      * 

This  system  is  identical  to  (6) ,  except  that   f   is  replaced  by 
f  cos(kd/2)   and  k   is  replaced  by  sin(kd/2)/ (d/2)  .   Thus  the 
solutions  for  u  (x,t)   and  h  (x,t)  ,  corresponding  to  Arakawa's 
initial  condition,  can  be  obtained  by  making  a  similar  replacement  in 
(12)  of  [5],  to  yield: 

2faVn   °°   .  .   .  v  sin(v  t) 

u  (x,t)  = /   Sin;aRJ —  cos(kx)  cos(kd/2)dk  ,      (13) 

c  it  Jn  ak  v 

0  c 

2HfaVn        °°   {l-cos(v    t)}       .     .  ..  •    ni/n\         n  a  in\ 

hc(t)»— jj-^     /     jS sinUk)    s.n(kx)    sxn(kd/2)cos(kd/2)    dk    ?       (u) 

c 

where 

v2     =   f2(l  +  A(X/d)2   sin2(kd/2))    .  (15) 

c 

Observe  that  (13)  is  not  obtainable  from  (9)  simply  by  replacing  v  by 

v   .   This  is  because  of  the   cos(kd/2)   term  that  arises  from  the 
c 

spatial  averaging  of  v  and  u   that  occurs  in  the  first  and  second 

equations  of  (11),  respectively.   Repeating  these  calculations  for  the 

other  schemes  will  show  that  u(x,t)   in  Scheme  D  also  requires  a 

cos(kd/2)   term,  and  that  h(x,t)   is  given,  in  general,  by 

2HfaV_   °°      ,-         ,   .  ,  .  . 
,  ,      .v        0   c  _,.  N  {1  -  cos  vt)  sm  ak    .  ,.     ...         .,,. 
h(x,t)  =  J   D(k)  7y r sm(kx)dk        (16) 


where 


2        ak 
0  v 


(6) 


D(k)  = 


and 


v/f  = 


sin(kd) 
d 

sin(kd/2) 
(d/2) 

cos(kd/2)sin(kd/2) 
(d/2) 

cos2(kd/2)sin(kd/2) 
(d/2) 

sin(kd//2) 

(d//2") 


2  2 

[1  +  aV] 


1/2 


[1  +  (A/d)2  sin2  (kd)] 


1/2 


[1  +  4(A/d)2  sin2  (kd/2)] 


1/2 


Differential 
Scheme  A 
Scheme  B 
Scheme  C 
Scheme  D 
Scheme  E 


,   Differential 
,   Scheme  A 
,   Scheme  B 


1/2 


[cos2  (kd/2)  +  4(A/d)2  sin2  (kd/2)]      ,   Scheme  C 

2  2    2      1/2 

[cos   (kd/2)  +  (A/d)   sin   (kd) ]        ,   Scheme  D 


[1  +  2(A/d)2  sin2  (kd//2~)] 


1/2 


,   Scheme  E 


(17) 


(18) 


In  Figure  2,  the  integrals  given  by  (16)  are  used  to  compute   (h/H)   at 
t  =  80  hours  in  the  same  manner  as  done  by  Winninghoff,  i.e.,  using 
Simpson's  rule  with  600  subintervals  in  the  interval   (0,Tr/d)  .   The 
relevant  parameter  values  used  were 

f  =  10~4  sec"1 


-2 


g  =  10  m  sec 


H  =  10  m 


a  =  d  =  A/2  . 


(7) 


These  graphs  too  should  be  compared  to  Figure  3,  which  shows  the 
equivalent  evaluation  of  (10),  with  only  v(k)   replaced  by  its  equiv- 
alent finite  difference  value.   (Figure  3  appears  identical  to  Figure 
3.6  in  Winninghoff . )   Note  especially  the  effect  on  the  solution  to 
Scheme  D.   We  conclude  that  Winninghoff  obtained  such  ill-behaved 
solutions  in  this  case  because  he  neglected  the  cos(kd/2)   term  in 
the  numerator  of  u(x,t)   in  these  cases,  and  this  term  helps  to 
control  the  singular  behavior  of   [v(k)]    as  k-*Tr   in  Scheme  D. 

Winninghoff  eliminated  Scheme  D  from  further  study  due  to  the 
extremely  noisy  behavior  shown  by  Figure  3.   Since  his  calculation  was 
erroneous,  it  can  not  be  discarded  for  this  reason,  however,  as  he 
showed,  it  is  still  seriously  deficient  in  terms  of  phase  speed  errors, 
and  should  probably  be  discarded  for  this  reason. 
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FOOTNOTE 

1.  A  sgn   term  was  erroneously  omitted  from  certain  expressions  in 

[5].   Thus,  in  [5],  equation  (13)  should  read, 

oo  I  i 

hg(x)  =  h(x,0)  -  -~-     !      sgn(x  -  s)e",X  "  s|A{|||(s,0)  -  v(s,0)}ds  , 
2  A  f  -oo 

the  equation  immediately  above  (28)  should  read 

00  —  8  x  A         o° 

1     f     iq     ikx      -  Ax  e   '  '  sgn(x)    v     ,     .      .   N 

2^   J        2  2  e    dk  =   ~T^2---—2Aj2      I        6(X  "  (2n  "  1)Ax)  ' 
1  +  A  a  A  [A  +  (Ax)  ]    n=-°° 


—  00 


equation  (30)  should  read 

H 


h  (X,0)  =  h(x,0)  +  9   1/9 

b  2Af[A  +  (Ax  )]  ' 


oo  I  i 

x  I        {sgn((2n  -  l)Ax)e"6'(2n  "  1>Axl/x  d(x  -  (2n  -  l)Ax,0) (2Ax) }  , 
n=-°o 

and  equation  (2.8)  should  be  identical  to  the  second  equation  in  this 

note. 


(9) 
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